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We discuss the measurement of new physics in long baseline neutrino oscillation experiments. 
Through the neutrino oscillation, the probability to detect the new physics effects such as flavor vio¬ 
lation is enhanced by the interference with the weak interaction. We carefully explain the situations 
that the interference can take place. Assuming a neutrino factory and an upgraded conventional 
beam, we estimate the feasibility to observe new physics numerically and point out that we can 
search new interactions using some channels, for example Vfj, ^ Vfj,, in these experiments. We also 
discuss several models which induce the effective interactions interfering with the weak interaction, 
and show that some new physics effects are large enough to be observed in the oscillation enhanced 
way. 

PACS numbers: 13.15.-|-g, 14.60.Pq, 14.60.St 


I. INTRODUCTION 


Neutrino oscillation induced by light neutrino masses and lepton mixings gives a plausible explanation for the results 
of the many neutrino experiments. In the three active neutrino framework, the neutrino mixing is given by, 

z/a = ^ UaiVi {a = e,pt,T), (1) 

i=l,2,3 

where a is the flavor index, i is the mass-eigenstate index. The mixing matrix [/ is a 3 x 3 unitary matrix, called 
Maki-Nakagawa-Sakata(MNS) matrix Q, defined as 

/ 1 0 0 \ / Ci3 0 \ / Ci2 Si2 0 \ 

[/ = I 0 C23 S23 I I 0 1 0 II —Sl2 Ci2 0 I , (2) 

\ 0 —S 23 C23 j Y —si 3 e~*^ 0 C23 y y 0 0 1 y 

where Sij{cij) stands for sin61y (cosdy). The atmospheric neutrino anomaly strongly suggests > Vt oscillation 
with large mixing, 4|C7^3p(l — |f7/i3p) > 0.9, and a larger mass squared difference, ~ 2.5 x 10“^eV^, which is 
almost confirmed by the K2K experiment and is expected to be reconfirmed in near future ^ . The solar neutrino 
deficit H] is also explained by the oscillation of Ve into another neutrino state. It gives several allowed regions for 

|CJei/t^e2| and the smaller mass squared difference among which the region for large mixing MSW 0 (LMSW) 

solution seems most preferable. To survey the LMSW region the KamLAND experiment will start soon. Also the 
Borexino experiment Q will strongly constrain the parameter region. 

However, we do not have sufficient information to determine the all mixing angles and the sign of the mass squared 
differences. |f/e 3 | is constrained strongly by CHOOZ Q and Palo Ver de jl^ experiments; |f7e3| < 0.16, but it has not 
determined yet. We also have no information about the sign of 5rn^^ |3l||, and the CP phase <5. Therefore there are 
many proposals for neutrino oscillation experiments to determine them. The future neutrino-oscillation experiments 
based on an accelerator, on both a conventional beam and a muon storage ring Q, are expected to give high 
precision tests of oscillation. Indeed according to these studies, we will determine the mixing angles and mass squared 
differences very precisely. The parameters relevant with atmospheric neutrino anomaly will be determined with error 
of a few %. In these experiments we can explore very small value (~ 0.01) of lie's a-nd we have a chance to observe 
the CP violation in the lepton sector. 

Till today main concern on future neutrino oscillation experiments is how precisely we can determine the oscillation 
parameters. Is it all we can do in such experiments ? The answer is NO.|^, In addition to the neutrino oscillation 
parameters we can probe the new lepton-flavor violating physics. It affects for example muon decay, matter effect 
and so on. Furthermore, as we will see in this paper, the effect of these exotic couplings is enhanced in oscillation 
phenomena. Even if the coupling is rather small, such new physics modifies the oscillation pattern distinguishably 
and hence it will be detectable. 
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In this paper we investigate the possibility to measure the effect of new physics in future neutrino oscillation 
experiments. We first study how the new physics contribute to the neutrino oscillation phenomena in section II. 
There we will see the reason why its effect is enhanced in the oscillation physics. In section III, we embody the 
formalism and adopt it for a neutrino factory and an upgraded conventional beam. We present numerical results 
on the feasibility to search new interactions in section IV. We argue that the energy dependence of event rate at 
relatively high energy region is a key issue for a observation. We will discuss the relation between lepton flavor 
violating interaction and some models in section V. Finally, a summary and discussion are given in section VI. 


II. BASIC IDEA 

In this section we reexamine what is measured in oscillation experiments and see how new physics gives contribution 
to the oscillation phenomena. Though we discuss only a neutrino factory to make the argument concrete, same 
discussion is followed in oscillation experiments with a conventional beam. 

First we remind ourselves what is really measured in a neutrino factory. All we know is that the muons, say, with 
negative charge decay at an accumulate ring and wrong sign muons are observed in a detector located at a length L 
away just after the time Ljc^ where c is the light speed. This is depicted schematically in Fig.|^. 



Since we know that there is the weak interaction process, we interpret such a wrong sign event as the evidence of 
the neutrino oscillation, Vf. —> z/^, which is graphically represented in Figj^. 



FIG. 2: Standard interpretation of a wrong sign event. 

Now if there is a flavor-changing exotic interaction, e.g., 

(3) 

then we will have the same signal of a wrong sign muon, whose diagram is shown in Fig.^, just like that caused by the 
weak interaction and the neutrino oscillation. We cannot distinguish these two kinds of contribution. The quantum 
mechanics tells us that in this case, to get a transition rate, we first sum up these amplitudes and then square the 
summation. Therefore there is an interference phenomenon between several amplitudes in this process. Through this 
interference we get an enhancement of the effect of new physics, that is, we can make oscillation-enhanced search for 
new interactions with neutrinos. 
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FIG. 3: Diagram which gives same signal as that given by Fig-i 


More formally we illustrate this situation as follows. Denote the transition amplitude between the initial state 
A + T and the final state C +U through an intermediate state B as where a final state U is one of 

the possible states which are unobserved like B. The transition probability between A + T and C +U, which is the 
absolute square of the sum of transition amplitudes with the same final state, is 


P{A + T ^C+U) 


2 


J2HAT;B-,C,U) 

B 


The transition probability from M to C is given by summing up all unobserved initial states T and final states hi as 


P[A^C)=Y^P{A + T ^C+U). 

r.u 


( 4 ) 


Furthermore suppose that the amplitude ^{A,Tq\Bq\C,Uq) is dominant over the amplitudes with intermediate states 
B{^ Bo) and/or initial states T/yf %) and/or final states hi{^ Uo). Then the dominant contribution to the probability 
P{A C)[^ is given by P{A + Tq ^ C + Uo) and deformed as 


P{A^C) ~ P{A + To^C+Uo) 
= \^{ATo-,Bo\CMo)? 

+ 2Re 


^ATo-,Bo-,CMe)* Y. ^{ATo-,B-CM^) 

2 

^ <t>iA,To-,B-,C,Uo) 

B^Bo 


( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 


Here the second term in eq.(|^) is given by the interference among the leading amplitude ^{A,To;Bo',C,Uo) and the 
sub leading amplitudes ^{A^To]B{^ Bo)]C,Uo)- Note that the interference arises among the processes which have 
quite the same initial and final states, even if some of them are unobserved. “Same state” means a state with not only 
the same particle species but also all other same physical quantities, such as same energy and same helicity. This is 
very important in the calculation of the transition rate. 

The leading amplitude is expected to arise from well known physics. On the other hand the amplitudes ^{A,To\B{^ 
Bo);C,Uo) is relatively small and may contain the contribution from new physics. However, even though \ ^{A,To\B{^ 
Bo)tC,Uo)\^ is very small, it would be possible for the interference term to be large. Thus even if the effect of new 
physics is not detectable due to systematic errors in a direct measurement, we may see the oscillation-enhanced effect 
of new physics. 

In neutrino oscillation experiments, the state A consists of the state of the parent particle of neutrinos and T consists 
from the state of the target particle that interacts with a neutrino in a detector. The final state C corresponds to the 
state of the particle which can induce an identifiable event in a detector as a result of the interaction with the target 
particle (T). Other final states U denote the states of all unobserved particles, which appear at both the neutrino- 
creation process and the detection process. The intermediate state B is interpreted as a neutrino state. The leading 
amplitude is given by the weak interaction and the neutrino oscillation, so the final state Uo is produced through the 
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weak interaction, and the intermediate state Bq is a nentrino state which “oscillates” from a certain weak interaction 
eigenstate to same/another one. Even if there is no oscillation, different intermediate states Bq) induced by 
exotic interactions give non-vanishing amplitudes for the total “oscillation” process. The new physics would violate 
the flavor conservation and/or the standard chiral property. Amplitudes which have all the same initial states Tq and 
final states lAa interfere with each other. It means that the leading amplitude from the weak interaction interferes 
with an exotic amplitude. In a neutrino-oscillation experiment such interference would give sub-leading contribution 
to the total probability. A large difference from the standard oscillation will be realized. 
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-I- . 

FIG. 5: The transition probability of a wrong-sign muon production process through the neutrino “oscillation”. We denote 
T and T' as the initial and the Hnal state of the target particle of the neutrino in a detector. Large dots describe effective 
interactions. Effective interactions characterized by Gf and Gt are induced by the weak interaction. We assume that the 
other effective couplings, like A®, are caused by new physics and their superscripts identify the types of effective interactions. 
The bold line in each diagram stands for the neutrino oscillation including all matter effects induced by the weak interaction 
and exotic interactions. Here, appropriate integrations of momenta and summations of helicity states for unobserved particles 
are omitted. 

To put the situation concretely, we consider a wrong-sign muon production process through the neutrino oscillation 
f'e —*■ In this case the initial state A is /i“ and the observed final state C is The above consideration tells that 

we have to calculate the transition rate depicted in Fig.^ Namely eq.(^ corresponds to Fig.^: The correspondences 
for the other symbols are T <-*■ T and U ^ D + T', respectively. In this process, the new physics effects contribute to 


















5 



2 


+ 2 Re 




Interference 0 ( 9 ) 

FIG. 6: The transition probability of a wrong-sign muon appearance process up to sub-leading contribution. 


three parts; the neutrino-creation process as mentioned above in eq. (i, the detection process, e.g., 

(9) 

and the matter effect described as 

^ee ^e/2 
^er 

where Ue is the electron number density in matter, For anti-neutrinos such an exotic matter effect is given by 
—H*. The interaction Hamiltonian of the exotic matter effect is induced in the similar way to the standard matter 
effect. 

By expanding Fig.^, the transition probability of this process is illustrated in Figj^ In Fig.^ the first diagram on 
the first line describes the contribution from the weak interaction and the oscillation process, where the unobserved 
state Ua corresponds to {e“, ,r'}. It gives the leading amplitude. The other diagrams on the first line in Fig.|| 

represent the effects of new physics, in which the states of neutrinos before and/or after the oscillation differ from the 
weak interaction case in flavor and/or chiral property. Since diagrams in the first line have completely the same final 
states C and Uq, they interfere with each other. The diagrams on the second and the third line in Fig|| have different 
final states form those of the first line, so that there is no interference with the weak interaction and it is expected 
that their contributions to the transition rate is negligibly small. 

Thus eq.(^) is given by the diagram of the first line in Fig.||. Figure || shows the eq. (i and eq. ©• The leading 
probability is proportional to since it is known that exotic effective interactions are so small according to 

the direct searches, that |Ap ^ \Gf\^ and \g\^ ^ \Gt\^■ However, the effect of new physics gives 0{X) and 0{g) 
contribution to the signal of wrong sign muon due to the interference with the amplitude of the weak interaction. 
This fact makes the search for new physics possible in an oscillation-enhanced way. It gives relatively large effect. On 
the contrary, even if the stored muon are very high intense, a direct detection of an exotic decay process will be very 
difficult due to systematics since the probabilities of such processes are proportional to the squares of the effective 
coupling constants and hence is expected too small to be detected. 


( 10 ) 













III. FORMALISM AND PARAMETERIZATION OF AN INTERFERENCE EFFECT WITH NEW 

PHYSICS IN NEUTRINO OSCILLATION 
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In this section we consider the neutrino oscillation in the presence of new physics more concretely to get the analytic 
expression for the transition rate. Hereafter we consider only diagrams which interfere with that of the known physics, 
le., the weak interaction. 

First we note that the amplitude for “neutrino oscillation” can be divided into three pieces: (1) Amplitude relevant 
with decay of a parent particle denoted as here C describes the type of interactions. For ^ decay, as we will see 
in eq. (0) and (^, there are two types of interactions, C = L,R while for tt decay we do not need this label, a 
distinguishes the particle species which easily propagates in the matter and make an interaction at a detector. (2) 
Amplitude representing a transition of these propagating particles, which are usually neutrinos, from one species a 
to another/same (3 , denoted as Tap. (3) Amplitude responsible for producing a charged lepton I from a propagated 
particle /3 at a detector, stood for by Here I denotes an interaction type. Using these notations we get the 

probability to observe a charged lepton at a detector as 


Pfj."^i+(i-) 


2 


Paf3D(3l± 

a/3CI 


a/SCI a’jS’C'I’ 


( 11 ) 


Therefore we can consider the effect of new physics separately for decay, propagation and detection processes. 

First we consider the decay process of parent particles. Since all final states must be the same, for a neutrino 
factory, the exotic decays of muons which are /i“ —> e~VaVe and /i“ —> can be amplified by the interference. 

Though in the presence of Majorana mass terms neutrinos and anti neutrinos can mix with each other, this effect 
is strongly suppressed by Therefore we do not have to consider decays into neutrino with opposite chirality 

such as fj,~ e~Vfj^Vj 3 The former is relevant with fi~ —> l~ and the latter is relevant with fj,~ —> l~^. In other words, 
we can approximate neutrinos to be massless except for the propagation process. This fact and Lorentz invariance 
allow only two kinds of new interactions in this process. For a wrong-sign mode, the allowed two interactions are the 
(U-A)(U-A)type, 

2^/2\a{v^YPL^J){e-fpPLVa)■, a = /T, T, (12) 

which has the same chiral property as the weak interaction but violates the flavor conservation, and the {V — A)(V+A) 
type, 

2V2X'a{i7pY^a){ejpPRn), a = e,/r,r. (13) 

The latter has different chiral property from the former, so that it gives different energy dependence to the transition 
rate. These exotic interactions interfere with the leading amplitude and contribute as next leading effects. Note that 
generally A and A' are complex numbers. |p^ 

In the case of the {V — A){V — A) type exotic interaction, we can introduce the interference effect by treating the 
initial state of oscillating neutrino as the superposition of all flavor eigenstates. On the pT —> pA process, we can 
take initial neutrino v as 


V = epVp+erVr, (14) 

where Cq, = Xa/Gp- This simple treatment is allowed only for the {V — A){V — A) type interaction because of the same 
interaction form as the weak interaction except for difference of the coupling constant and the flavor of antineutrino. 
In this case we can generalize the initial neutrino for any flavor, using Y. Grossman’s source state notation 0 , as, 

0 


= a,P = e,n,T, 


W = 


1 e'* e® 

^e/i. ^er 


We can include the total exotic effect into the oscillation probability as 


P. 








(15) 


(16) 
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This treatment is also valid for the effect on the oscillation. 

In the case of the {V—A){V+A) type exotic interaction, we cannot treat interference terms simply. The interference 
term between the weak interaction and an exotic interaction eq. (0) denoted as P^}^i gives the rate for the observation 
of the wrong sign charged lepton, which is interpreted normally as the oscillation from Vg —> vi in decay, as follows: 


P 


( 1 ) _ 1 + 1 r <Ppe _ ) 

-->1+ ~ o 2^ ! r, TP r, TP ° \PPv Pe Pv^j 


2 (27r)2 


spin 


X 2Re 




01 


ol(3' I' 


1 + P.8Gf 


a/3/3'//' 


(17) 


where Pp, is the polarization of the initial //“, E^ is v energy, and Pp(Ep) is // momentum (energy). 

In the case for the observation of the same sign charged lepton, which is interpreted as the oscillation from i/p —> i// 
in decay, as follows: 


P 


( 1 ) 


1 - Pa SGf 


+/- 


mempEai\Pp\ - Ep) 


a/3/3'//' 


(18) 


For TT decay the situation is much simpler. In the presence of new physics there may be a flavor violating decay of 
TT such as 7r“ ^ ^~Ua{a = e,T). This effect changes the initial v state; 

^tr+^pr^T- (19) 

In this case we do not have to worry about the type of new physics which gives a flavor changing tt decay at a low 
energy scale. Due to kinematics, the energy and the helicity of the decaying particles, p and v are fixed. 

Next we consider the propagation process. Exotic interactions also modify the Hamiltonian for neutrino propagation 
as M, 


Ppa — 


1 


^0 


17/3, 


Sml^ 


U, 




/a + Qee 


^er \ 1 


+ aep* 

cipp 

1 / ; 

(20) 

\ aer* 

* 

/3a J 


Uaf} is the 

extra 

matter effect 

due to new physics 


interactions, that is defined by Uap = 2^/2t'^pGFUeE^■ Note that to consider the magnitude of the matter effect, the 
type of the interaction is irrelevant since in matter particles are at rest and hence the dependence on the chirality is 
averaged out.fl^ 

Finally we make a comment about new physics which affect a detection process. To consider this process we need 
the similar treatment to that at the decay process, that is, we have to separate contribution of new interactions 
following the difference of the chirality dependence. However to take into account new physics at a detector, the 
parton distribution and a knowledge about hadronization are necessary. Though we may wonder whether we can 
parameterize the effect of new physics at the detector g/Gr as like e®. It is expected that has a complicated 
energy dependence due to the parton distribution for example in a energy region of a neutrino factory. Consider 
the case that there is an elementary process from lepton flavor violating new physics including strange quark. To 
parameterize its effect we need both its magnitude and the distribution function of strange quark in nucleon which 
will show the dependence on the neutrino energy (more exactly, the transfered momentum from neutrino to strange 
quark). They are beyond our ability and hence we do not consider them further in this paper, though new physics 
which can affect the decay process have contribution to the detection process too. 


IV. NUMERICAL ANALYSIS 

Let us discuss the feasibility to observe the signal induced by new physics. We will deal with both an upgraded 
conventional beam and a neutrino factory but present a little bit qualitatively different analysis in each experiment. 








For a neutrino factory, we use following procedures: First, we assume the magnitude of the effects caused by new 
interactions. Next, we calculate the event numbers including the effects of new physics and also calculate that 

based on the standard model . Then, we define the following quantity, so called function |Q, 


bin 


X' 


-E 






( 21 ) 


where i is the energy bin index, is the muon number, and M(jet is the detector mass. To claim that new physics 
effects can be observed at 90% confidence level, it is required 


X^ > X90%> 


and this condition is rewritten as 




X90% 

Y2 

fact 


( 22 ) 


(23) 


From the above method, we can obtain the necessary muon number and the detector mass to observe the new physics 
effects at 90% confidence level. 

On the other hand, there are some kinds of option on beam configurations for an upgraded conventional beam; 
wide band, narrow band, off axis, and so on. We discuss the event number of no-oscillated neutrino events at the 
detector Ni, unlike the case for a neutrino factory. The concrete procedure is almost the same as that for a neutrino 
factory. We separate function defined in first line of eq.(^) into two parts; 


X^ = 


(24) 


and get the necessary number of no-oscillated neutrino events from 


N, > 


X90% 

V2 

-^conv 


(25) 


In the numerical calculation for a conventional beam, we consider the wide-band-beam like situation, whose flux 
distribution for energy is constant. 


A. (V — A){V — A) type new interaction 


Here, we deal with the case that there are only (V — A)(V — A) type new interactions in the lepton sector. In 
this case, we need to consider the effect represented in eqs.(p^ and (20). Before making the presentation of the 
numerical calculations, we give the analytic expression for the sensitivities to understand the essential features. As 
we showed in section III, the interference terms between (V — A){V — A) type interactions and the weak interaction 
have the same dependence on the ^ polarization leading term does. It can not be expected that the sensitivity to such 
interference term becomes better by the control of the parent ^ polarization. We consider here an unpolarized muon 
beam. The “oscillation probability” is given by eq. (0) in this situation. More detailed calculations are presented in 
the Appendix. 


1. r'e —> channel in a neutrino factory 

The analytic expression of probability for v,, given in the Appendix A shows that the effect due to and 

e®’™ are irrelevant since these terms are proportional to sin^ 20 i 3 x e in the high energy region, so it is difficult to 
observe their effects. The flavor changing processes between muon and electron, e.g., fi ey, /r e conversion, are 
strictly constrained from experiments, and as we argue in the next section the box diagrams of the ^-to-e processes 
must relate to and e®^. Therefore, the magnitude of e®^ has very severe bound and the terms depending on it are 
also not effective. Assuming some models, e.g., MSSM with right-handed neutrinos, it is expected that the magnitude 
^a/3 ^ 0/3 same order because they are produced by similar diagrams. This is also discussed in the next 

section. We assume naively that the terms depending on are constrained as e®^. On the other hand, the r-to-e 
processes do not give the tight bound to Hence, we investigate the effect induced by e®^™ hrst. 
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Before surveying the required N^M^et for each baseline L and muon energy £1^, we see the behavior of contribution 
of to the “oscillation probability” to consider the optimum setup for L and In the high energy region such 
as the matter effect a is much greater than Sin'll, the first order contribution of to the transition probability, 
is constructed by four parts that have the different e®^™ dependences: 


— 2s23S2x23S2x13 

„2 


a , \ / Sm^ 


31 


Ci3 (s5Re[e®^] - C5lm[e®^]) 


+ c?3 (c5Re[e®^] + S5lm[e®^]) 1 


- c ?3 (s,Re[e™] + c,Im[e™J) f^L 



- o«i3 (c5R.e[e™] - S5lm[e™]) 


4£;. 


L +2 


dm? 


31 


4£;. 


4£;. 


(26a) 


L (26b) 


Sm 


31 


4£;. 


(26c) 

(26d) 


where S 2 xij = sin20y. 


Since we can suppose that is proportional to Ef^ in a neutrino factory, E^ and L dependence 


of the sensitivity to each term can be approximated as 


X^(26b)cx <^1-- 



2 V4-E, 


X E,, 




oc 


4£;„ 


L] xE^, x"(26d)cx 


iE„ 


L X E„ 


(27) 


All of eq.lCT) are proportional to £1^, so the sensitivities must get better as the energy becomes higher. 


Each of eq.(27) depends on L in different way. (p6b|) become tiny for longer baseline length within the region 


that we are now interested in. This fact means that a shorter baseline experiment has an advantage over a longer one 
to obser ve (|26bD ’s ef fect. By contrast with this, it is found that longer baseline will be better to search for the effects 
of (^), ( p!6cj ) and (26d). 

Each of eq.(p^ depends on the combinations of e®^™ and the CP phase <5. What we observe is the combination of 
them. The effects of e’s can be sources of the CP-violation effect. ||l7|] In the discussion about the observation of the 
CP phase, this fact should be considered. The analytic expressions also show that the sensitivities are proportional 
to |e|2. 

Now, we show the results of the numerical calculations. The parameters that we use here are 


smfyi2 = 


1 

2 ’ 


sin023 = sindi3 = 0.1, 

v2 


Smh = 5 X 10 


5mh = 3 X 10 

TT 

2 ’ 


(28) 


and take |e| = 3 x 10“^, which is a reference value for the feasibility to observe the effect by using the method of 
the oscillation enhancement. Except for e®^ and e®g, the constraints of the processes of charged lepton have not 
forbidden this magnitude of e’s. 

Fig.^ shows the required TV^Mdet in the case where we do not take into account the uncertainties of the mixing 
parameters. We can check whether the approximated equations, eq.(|^), are correct from the behavior of the plots. 
As eq.(p^ implies, contribution from the new interaction depends on the combinations of e and S. We take S = 7r/2, 
so we can extract each term of eq.(|^) by taking e®^™ pure real and im agina ry. Th erefo re, the plots of Figj^from left 
to right correspond to the required data size to observe each term in eq.(26a) to e q.(^6d| ) respectively. The behavior of 
these plots is consistent with the expectations from the analytic expressions in eq.(P7D. To consider realistic situations, 
the new physics effects in both source and matter must be taken into account simultaneously. We present Fig|| as 
the plots in the cases where the same magnitude of new effects in both source and matter are introduced. They show 
that the total effects are given by the simple summation of each effect. 

In realistic situations the uncertainties of the theoretical parameters have to be taken into account.Once the 


uncertainties are introduced, it can be expected that the sensitivities shown in Figs.^ and 
The e®^™ effect can be absorbed easily into the main (unperturbed) part of oscillation 


will be spoiled completely. 
3y adjusting the theoretical 


parameters since the effects have the same energy dependence as the main part has. Indeed, taking into account these 
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S^[GeV] 




£^M[GeV] 


^^^[GeV] 


FIG. 7: Contour plots of the required A^^Mdet to observe the new physics effects concerning at 90% C.L. in ^ r';/ channel 
using a neutrino factory. From left to right: (eerjCf^) = (3.0 x 10“®,0), (3.0 x 10“®i,0), (0,3.0 x 10“®), (0,3.0 x 10“®i). The 
uncertainties of theoretical parameters are not considered in these plots. As we point out in the text, when the uncertainties are 
taken into account, the sensitivities are completely lost. In Fig. to Figjlll contours mean 0.01, 0.05, 0.1, 0.5, 1x10^^ • lOOkt. 



F:^[GeV] 


SmIGcV] 


SmIGcV] 


F^^iGeV] 


FIG. 8 : Same as Fig.M, but here new interactions in both source and matter are taken into account. From left to right and top 

, o o. , o o., o o . . o o. 


to bottom: (eer,«(((-) = (3.0 x 10 ,3.0 x 10 ), 


(3.0 X 10“®i,3.0 X 10“®ij 


(3.0 X 10“®, -3.0 X 10“®), (3.0 X 10“®i, -3.0 x 10“®i), (3.0 x 10“ 
be understood by addition or subtraction of the plots in Fig.|^ 


(3.0 x 10“®,3.0 X 10“®i), (3.0 x 10“®i, 


3.0 X 10“®); 


, —3.0 X 10 ®i), (3.0 X 10 ®i, —3.0 x 10 ®). These plots can 


uncertainties, the sensitivities to are completely washed out. Therefore, we have to look for the terms whose 
energy dependence differ from that of the main oscillation term in the high energy region. In i/g ^ channel, the 
effects caused by have such energy dependence. It can be represented analytically in the high energy region as 


— 2s23S2x13 

(s5Re[e®^] - C5lin[e®^]) 


X |l 3 ( 4 ^/) + 3 (2C2X13 H 3 CI 3 ) ( 4 ^/ 
- (c5R-e[eep] + S5lm[e®^]) 


Sr^ 


6r^ 
4E, ^ 


(29a) 
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+ 2c^ 3 (s5Re[e™ ] + Qlm[e™ ]) ( -^L 
+ 2 (c5Re[e™] - S5lm[e™]) 


' Smi 


31 


^ r ■ ^ +^4^2x13) 


(29b) 

(29c) 



(29d) 


Contribution for the transition probability labeled (29a), ( |29b| ) and ( P9d| ) depends on IjE^. Consequently, the 
sensitivities to the terms must be robust against the uncertainties of the theoretical parameters since in the high 
energy region they can be distinguished from the main oscillation part by observing the energy dependence. The 
claims mentioned above are confirmed numerically by Fig. i| an d |l0| . By comparison of these graphs, we can see 
that the sensitivities to observe the contribution of (29a), (291) and (29d) do not suffer from the uncertainties. |3^ 
Incidentally, we note that though the uncertainties wreck the sensitivity to (29c) since it is proportional to 1/4 > 


second order term brings constant contribution for energy and this signal does not vanish. 



FIG. 9: Contour plots of the required A^^^Mdet to observe the new physics effects concerning at 90% C.L. in —> Vfj, 


channel when there is no uncertainty for theoretical parameters. From left to righ t: fe f,,., 
(0,3.0 X 10”^), (0,3.0 X 10“^i). Each plot corresponds to the sensitivities to eq.(29a)~eq.(b 


^ (3.0 X 10-®, 0), (3.0 X lO-^i, 0), 
respectively. 



FIG. 10: Same as 


Fig.H, but here each parameter has 10% uncertainty. 


2 . Vfj, —> Up, channel in a neutrino factory 

Only the effect that depend on e®’™ will be large enough to be observed in Vp disappearance channel. As we 

show in the next section, this quantity is not strongly bound by the charged lepton processes. The analytic expressions 
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for the terms concerning e®’™ (see the Appendix B) indicate two facts: (1) this channel sensitive only to the real part of 
and (2) the effect that comes from the real part of will be small in the assumed parameter region. In addition, 
it shows that the terms depending on e®’™ are hard to be absorbed by the uncertainty of the theoretical parameters. 
Note that these terms do not depend on the CP phase 5, that is, it is expected that we can get information on the 
phases of e^’™. The sensitivity plots calculated numerically are shown in Figji]. They behave as expected by the 
analytic expressions. The uncertainties of the theoretical parameters do not affect the observability since the terms 
that depend on do not vanish in the high energy region. The sensitivities depend strongly on the phase of e®’™ 

but do not depend on 5. We can directly know the phase of the lepton-flavor violation process without care of 5. 



FIG. 11: Contour plots of the required A^Mdet to observe the new physics effects concerning e®’™ at 90% C.L. in —> Vfj, 
channel. All theoretical parameters are assumed to have 10% uncertainty. From left to right, each plot corresponds to 
«r,c;rr) = (3.0 X 10-^0), (3.0 X 10-®i,0), (0,3.0 x lO"®), (0,3.0 x 10“®!). 


3. Channels with r and e observation in a neutrino factory 

The technologies for tau observation in detection and the charge identification of electron to distinguish i/^ 
with Uf. 1^^ are under R & D. If it is possible to observe these particles clearly, what can we get ? 

In Ve Vt channel, we can explore e®;!" (see Fig.^^. The uncertainties of the theoretical parameters will not 
disturb the observability. In ^ Vr channel, all we can observe is only the effect of e®’™ (see Fig.^^. In comparison 
with channel, we will not have so much benefit in terms of sensitivities to the magnitude of e.|^^ However 

in this channel unlike channel, the observable must depend on the combination of the e®’™’s phase and the 

CP phase d, that is, the observation in these two channels are qualitatively different. In z^^ —> channel, we can 
search the effect of not e®^ but e® g at the muon decay process and the effect of e™ at the propagation process just 
like z^e channel (see Fig. 0- In Vg —> i/g disappearance channel, oscillation effects themselves are much smaller 

than the no-oscillation signal. Though some effects of new physics give the different energy dependence from the main 
oscillation term, we will not be able to get any information for oscillation-enhanced new physics. 



FIG. 12: Contour plots of the required A^Mdet to observe the new physics effects concerning e®;-"* at 90% C.L. in —> Vt 
channel. All theoretical parameters are assumed to have 10% uncertainty. From left to right, each plot corresponds to 
{t%^,eT-r) = (3.0 X 10■^0), (3.0 x 10"^i,0), (0,3.0 x 10“®), (0,3.0 x 10“®i). 
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S^[GeV] 


£^M[GeV] 


£^M[GeV] 


E^[GeV] 


FIG. 13: Same as Figj^, but —> Ur channel. 



FIG. 14: Contour plots of the required A^/^Mdet to observe the new physics effects concerning and at 90% C.L. in 
channel. All theoretical parameters are assumed to have 10% uncertainty. From left to right, each plot corresponds 
to = (3.0 X 10"^0), (3.0 X 10"®i,0), (0,3.0 x 10"^), (0,3.0 x 10"®i). 


B. {y — A){V + A) type interactions 

The signals of {V — A) {V + A) type interactions are hard to be observed because they are suppressed by the factor, 
The interference term between the leading term and (V — A)(V + A) type ones has different polarization 
dependence from that of the leading contribution unlike {V — A) (V — A) type new interaction. If we can make goo d 
use of this fact, then we may be able to expect to gain somewhat better sensitivity. As we saw in eqs.(|^ and (|l j), 
the utilization of the muon polarization works only in Va- Figure jl^ describes the sensitivity to e'X'/G f) 

in Vfj, —> Vr channel. We can gain a little, but have no advantage over a direct search against our expectation. 



FIG. 15: Gontour plots of the required A^Mdet to observe the {V — A) {V + A) type new interaction, = 3.0 x 10“®i, at 90% 
C.L. in F/j, —> Ur channel. Gontours mean 1, 5, 10, 50, 100x10^^ ■ lOOkt. Here, the uncertainties of the theoretical parameters 
are not considered. From left to right, each plot corresponds to =0.99, 0, —0.99. 
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C. 


channel in an upgraded conventional beam 


In an experiment using a conventional beam, there are some different points from the case of a neutrino factory. 
Since we consider two body decay of tt’s, the created charged lepton and the neutrino have the fixed helicity whatever 
new physics is. It means that we do not have to worry the type of a new interaction. We can always parameterize the 
effect of new physics concerning tt decay by and When we consider a new source interaction, new interactions 
between neutrino and nucleon must be introduced not only in matter but also in detector. However to consider such 
effects we have to treat complicated hadronic processes, so we do not consider the relation between these two effects. 
The contribution to the source and the matter effect is not always generated by the same new interactions. Hence, 
we can suppose that different magnitude interaction in each process. 

As it became clear in analysis for a neutrino factory, all we can probe is the contribution of e®^, e™ in 
channel and that of e®’™ in channel. However, new interactions in matter can not be observed in the energy 


and the baseline regions that we assume here and hence we do not study their effect. The numerical results for the 
necessary number of no-oscillated neutrinos Ni, are shown in Figj^, Fig, 17, and Fig. E The neutrino event number 
of current proposed experiments is estimated of 0(10^ ~ 10^), and C>(1(F) times larger exposure is expected in the 
next generation. 

Figures and E show that the sensitivity strongly depends on the complex phase of e’s. This fact means that 

the observations give us information of the phase in lepton-flavor-violation interactions. Depending on models, some 
interesting issues are revealed. In the next section, we discuss models which may give the significant e’s. As we point 
out above, e’s for a neutrino factory and a conventional beam have different dependence on new interactions. Hence, 
it is important for new physics search to compare the results of a neutrino factory and that of a conventional beam. 



F^GeV] 


F^^GeV] 


FIG. 16: Gontour plots of the required no-oscillated neutrino number, A„, to observe the new physics effects concerning 
at 90% C.L. in ^ channel using an upgraded conventional beam. All theoretical parameters are assumed to have 10% 
uncertainty. From left to right, each plot corresponds to = 3.0 x 10“®, 3.0 x 10“^i. In this energy and baseline region, a 
experiment does not sensitive to new interactions in matter. In Fig. lil, Fig. 0 and Fig.[l^, contours mean 2, 4, 6, 8, 10 xlO^ 
no-oscillated neutrinos. 
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FIG. 17: Contour plots of the required no-oscillated neutrino number, to observe the new physics effects concerning 
at 90% C.L. in channel using an upgraded conventional beam. All theoretical parameters are assumed to have 10% 

uncertainty. From left to right, each plot corresponds to = 3.0 x 10“^, 3.0 x 10“^L 
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FIG. 18: Same as Figjl^, but —> v.^ channel. 

D. short summary for numerical analysis 

For a neutrino factory: 

• In t'a —!■ 12 ^, (a = e, /3 = e, /r, r) appearance channel, the observable effects of new {V — A){V — A) interactions 
come only from e*’^. The others are too small or too vulnerable against the adjustment of the theoretical 
parameters. Note that 5 and e’s phase are correlated. Namely the measured values are a certain combination 
of i5 and e. 

• In —!■ disappearance channel, we can measure e®’™ depending on their phase. In other words, the signal 

includes information of the phase. Furthermore, there is no correlation between S and e, so the measurement 
tells us directly the phase concerning the lepton-flavor violating process. In i/g —> disappearance channel, we 

can not get anything for new interactions in the oscillation enhanced way. 

• The is proportional to |ep. The expected sensitivity is to |e| > 0(10“"^) by using this methodology. 

• When the situations that new interactions exist not only in the source but also in the matter effect are considered, 
we can easily understand the sensitivity by simply adding each effect. 

• Oscillation-enhanced effects for the (V — A)(V + A) type interactions are strongly suppressed by We/m^, so we 
can not get an advantage over a direct measurement. 

For an upgraded conventional beam: 

• We do not have to care the types of new interactions in the source. The analyses for the feasibility are similar 
to that of {V — A){V — A) type for a neutrino factory. In the assumed energy and baseline region, there is no 
sensitivity to the new effect in matter. 

• The e’s for a conventional beam have different dependence from those for a neutrino factory on new interactions. 
Therefore, the comparison between two methods makes clear the species of new physics. 



V. INDICATIONS TO EXOTIC DECAYS FROM VARIOUS MODELS 

In this section we discuss models which give exotic interactions interfering with the weak interaction in the neutrino 
oscillation and survey to which process those exotic interactions contribute. We should consider models which give 
an explanation for the smallness of the neutrino masses and the lepton mixings. In those models we can expect that 
flavor violating processes are induced, and since to explain the neutrino masses and the lepton mixings we need to 
introduce sources of flavor violations. 

There are well-known two types of models which explain the smallness of the neutrino masses and the lepton 
mixings. One is a seesaw type and the other is a radiative type |2^. These models induce other flavor violating 
effects generally. Such other flavor violating contribution could be large enough to give observable effects for neutrino 
oscillation experiments through the interfering effects. We now consider these two types of models concerning mainly 
on supersymmetric models and study how large the exotic effective interactions are. 





16 



eVr, (4e) 
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bound from charged lepton processes 
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X 

O 

A' 

lAfi > Vt 

X 

A^ 

0 “ 



X 

X 

i/g ^ Ue 

X 

X 

X 


“appropriate mode but bound is too strict 

^too vulnerable against the adjustment of the theoretical parameters 
“depending on the e’s phase, in other words, sensitive to the phase 


TABLE I: Summary on the feasibility of (V — A){V — A) type new interactions. 


A. Supersymmetric models with right-handed neutrinos 


Among the promising extensions of the standard model, a supersymmetric standard (or grand unified) model with 
right-handed neutrinos is often considered. In this class of models, if the gravity-induced supersymmetry breaking is 
employed, through the renormalization effect, the large flavor violating slepton masses are induced |^, even if at the 
cutoff scale there is no flavor violating slepton mass. 



( 1 ) ( 2 ) 

FIG. 19: (1) An Example of one-loop diagrams contributing to the decay and the matter effect. (2) An Example of tree 
diagrams contributing to the tt decay, the matter effect and the detection process. Where Srrij denotes an off-diagonal element 
of the mass matrix for sleptons which violates the lepton flavor. 

Due to the lepton flavor violating masses, one-loop box diagrams constructed by propagations of superpartners 
induce exotic effective (V — A)(V — A) interactions, which contribute to the decay side, the matter effect and the 
detection side all together. One of effective interactions contributing to the decay side in oscillation experiments based 
on a muon beam is drawn in Figp^(l). This gives e®,.. There is a constraint on the magnitude of Figp^(l) from the 
lepton flavor violating decays of r such as r ^ /ry and r ^ jiee. However, these constraints are not so strict and 
can be of 0(10“^). This type of the effective interaction also induces the matter effect similarly replacing Ve with e 
and /i with respectively, and also can be of 0(10“^). By changing the external legs in Fig [T^(l), we can draw 
similar diagrams for four-lepton couplings which violate the lepton flavor, and using them we can estimate how large 
values the are in each model and experimentally. They can contribute to the fjL decay and the matter effects. 

In addition, there are effective interactions generated by one-loop diagrams including squark propagators as Fi^ 
(2). These kinds of diagrams can contribute to the detection processes and the matter effects. Furthermore they can 
contribute to the tt decay. 

These two kinds of diagrams affect the oscillation in a neutrino factory and a conventional beam differently. By 
comparing the results of these two experiments, we may have some information on the scalar masses. For example, if 
the squark masses are much heavier than the slepton masses, then only the figures of the type Figp^(l) contribute to 
the oscillation phenomena, and hence there is a difference between “oscillation probabilities” in these two experiments. 
On the other hand if the magnitudes of the scalar masses are comparable, these two experiments are affected from 
new physics in the decay side, the matter effect and the detection side similarly. 
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It is important that we can get information on the phases of e’s. Their phases arise from the phases of the flavor 
changing slepton masses, whose origins are the Yukawa couplings for vi^ and Vfi at the high energy scale. We may be 
able to acquire information on the phases of Dirac-Yukawa couplings. 


B. Models with radiatively induced neutrino masses 




Ir,l i 

I 

I 

-i-- 

e 


d! 


dR 


d 


( 1 ) 


( 2 ) 


FIG. 20: (1) Examples of tree diagrams contributing to the p decay and the matter effect. (2) An Example of tree diagrams 
contribnting to the tt decay, the matter effect and the detection process. The first diagram (1) induces both {V — A){V — A) type 
and {V — A){V + A) type interactions for the muon decay and the matter effect. The second diagram (2) induces {V — A){V + A) 
type interactions only, which contribute to the matter effect. 

There are many kinds of models in which the neutrino masses and the lepton mixings are radiatively induced. 
Among them we consider supersymmetric models with R-parity violating terms. A most general superpotential which 
breaks R-parity is 


IArpv — ^ajS-fLocL/sEy + X(y.jj^LctQi^D'^ A- -I- (30) 

where = ~A^o, 7 , A"^..^ = —A"..^^. La, E^, Qa, Ua, and Hua are the superfields corresponding to the 

lepton doublets, the right-handed charged leptons, the quark doublets, the right-handed up-type-quarks, the down- 
type quarks and the up-type Higgs doublet respectively. Three trilinear terms are allowed if the R-parity is broken 
generally. However there is no reason why all R-parity violating terms should be included at the same time. For 
example, we forbid interactions proportional to A”^^ by imposing the baryon-number conservation to avoid rapid 
proton decay. To induce small masses for neutrinos, the first term or the second term in eq.(|^) must be included. 
For example the first term in eq.(|^) has interactions described as 

£ = XaAjii'aL^lSLet^R + e^Rl^aLef3L + e^Re/SL^aL - {a ^ (3)) -I- H.c. (31) 

These interactions work in a same way as those of the Zee model. Namely, sleptons interact with charged leptons 
and neutrinos as the Zee boson, and the small neutrino masses and the lepton mixings are induced radiatively. Such 
interactions induce {V —^{V — A) type and {V — A) {V + A) type effective four-lepton interactions which violate 
the lepton flavor (see FigJM). These affect oscillation experiments through the interference with the weak interaction. 
Since there are only four-lepton type effective interactions, we would be able to see the difference between “oscillation 
probabilities” in two different experiments. In the case of a neutrino factory, they affect both the decay side and the 
matter effect all together, while in the case of a conventional beam they affect only the matter effect. From these 
kinds of Fig.|2^(l) e*’™ and are induced. In this case the magnitudes of e’s could be larger than the case in the 
previous subsection. 

Interactions generated by the second term in eq.(^), are described as 

“ ^a/3'f^l3S^i^aRdsLdjR -f d^RVahdsL + d^RdsL^ah) 

A ^a^'^iRpL^aLd’^R + djReaLUjSL + dr^RUpLeah) + H.C., (32) 

where is the mixing matrix for the quark sector. In this case squarks work as the Zee boson similarly. These 

interactions induce effective four-Fermi interactions with the lepton flavor violation, as drawn in Fig.|20|-(2). They 
affect the matter effect and the detection process in both experiments based on a muon beam and on a conventional 
beam. For a conventional beam, there is an additional interference in the decay side. Therefore if the R-parity 
violating term is limited only to the second term in eq.(^0|), then we also expect to be able to make sure the existence 
of such terms by comparison with the two different experiments. 
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C. General properties of the effective four-lepton interactions and order estimation for their couplings. 


Without assuming physics that gives the effective four-lepton interactions containing two left-handed neutrinos 
and two charged leptons with the lepton flavor violation, we can classify two types of such interactions from Lorentz 
invariance as the {V — A){y — A) type and the {V — A) [V -\- A) type generally. 

The flavor violating four-lepton couplings with the {V — A){V — A) type are classified into two categories: SU(2)j;, 
singlet type, 

hal3-/5iiaClf3){l^CHs) (33) 

= -j^haf3~fshaf3~i5{^aLCvi3L — f'aLC'e/3i)(e^iC'l'r'5L — 

and triplet type, 

9c.0A^c.T^Ci0){l^C^T-ls) (34) 

= gajS'iS {{eaLChipL + VahCej3L){e-fLC^VgL + ^Sl) 

-\-2{VaLCvi3L){y-lLC^VS l) + ‘^{eaLCei3L){e^LC^G sl)} ■ 

Here la is a lepton doublet with flavor a and r“ is a Pauli matrix. Due to the SU(2)i, invariance coupling constants 
must satisfy hap-yS = —hpajS = —haps-y and Qap-yS = gpa~fS = gpo-Sj- Above two interactions are described as the 
effective interactions induced by exchange of not only vector particles but also scalar particles. The former type is 
induced, for example, if there is a coupling of the form 

ICl4>s+il.c., (35) 

where 0s is an SU(2 )l singlet. The latter is induced by exchange of an SU(2)j;, triplet (j)^, if the coupling 

H.C., (36) 


exists. 

On the other hand, the flavor violating four-lepton couplings with {V — A) {V + A) type have only one category: 
SU(2 )l doublet type, 

fap^5{loi^PR){p^Rls') (iii^) 

= faP'fS {{eaLepR){e^ResL) + {yaLepR){e^RvsL)} ■ 

This type is induced by the exchange of a scalar particle belonging to an SU(2)i, doublet.!^ Suppose that there is a 
coupling of the form 


Ze0D + H.C., 


(38) 


where 0 d is an SU(2 )l doublet scalar, the interaction eq.(^) is induced. 

In general the constraint on the magnitude of the singlet type effecti ve co upling hap-yS is rather weak. A radiative 
model often contains such a singlet scalar. Indeed models in section VB contain such SU(2)i singlet scalars (see 
Fig.|^(l) with Ir). As for the p, decays which contribute to the interference phenomena, replacing a, 7 and S in 
hap~fS with e, p, e and r we get and replacing a, (3, 7 and <5 in hap^s with e, p, p and r we get e®,.. In this case 
there is no constraint from r — > III decays. Also constraints from t ^ p'y (for the former) and r — > ey (for the latter) 
are not severe. |2^ There are constraints from the universality and it gives the stringent constraints on them. The 
magnitudes of e’s can be 0(10"^) using the universality which is employed in Ref, . However this constraint 
can be relax ed significantly. |^, Thus we can expect rather large e®’™ in this kind of models. That is, models in 
section |VB| would give very large e®’™’s. 

For the triplet type and the doublet type there are stronger constraints. The models in section VB realize the 
doublet type (see Fig.p0[-(1) with Ir). Supersymmetric models with right-handed neutrinos in sectioiT~ |V^ gi ve both 
the triplet type and the singlet type interactions. In these cases the discussion by the authors in Ref.|[17|| can be 
applied ETl, that is, 


e®^ < 3.1 X I0“^ e®^ < 5 X I0“®, e®^ < 3.2 x 10”^ 


(39) 
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VI. SUMMARY 

In this paper we considered how we can observe the effect of new physics in a neutrino oscillation experiment. 

First we reminded ourselves what we really measure in a neutrino oscillation experiment. All we really observe 
are a decay of parent particle (tt or /r) at an accelerator cite and an appearance of a charged lepton at a detector. 
Neutrinos behave as merely intermediate states. Therefore if new physics gives the same amplitudes as those given 
by the weak interaction, we cannot distinguish those amplitudes and hence the transition amplitudes given by new 
physics can interfere with those of the weak interaction. It means the effect of new physics can be amplified quantum- 
mechanically. That is, we have the effect of not 0{X^) but 0(A), where A is an effective coupling of new physics, in 
an oscillation experiment. We understood also that the all particle states which appear as external lines including 
unobserved particles must be the same for the interference to occur. This statement means that not only the particle 
species but also their energy and helicity must be the same. We have to take into account both the particle species 
and a type of an interaction. 

Next we derived the transition rate for an appearance of a charged lepton, which is usually interpreted by the 
neutrino oscillation. To calculate the transition rate in a neutrino factory, we have to consider not only particle states 
in external lines but also a type of new physics. We separated the effects of new physics in three parts: (I) We 
considered new physics affecting the decay process of a parent particle. If new physics which affects the decay is a 
{V - A) {V- A) ty^, eq.Q, then it changes the initial state of neutrinos from a pure flavor eigenstate to a mixed 
state given by eq.(^. On the contrary, for the case of a (U — A)(y + A) type new interaction, eq.(p^), we have a 
rather complicated equation for the transition rate, though this effect is suppressed by me/m^ due to the difference 
of the chirality dependence on the interaction from that on the weak interaction. On the other hand, in an oscillation 
experiment with a conventional beam, we do not have to worry about the interaction type for the tt decay, since 
the energy and the helicity of mouns and neutrinos are fixed due to the kinematics of the two body decay. We can 
parameterize the effect of new physics using eq. dil)- (2) We studied the effects of new physics in matter. They are 
parametrized as eq. ( 0 ) for both oscillation experiments. Their effects also appear linearly in and hence they 
can give significant modihcation on an “oscillation event”. (3) We gave a comment how new physics appears at the 
detector. The essence is quite similar with that of the decay process of a parent particle. However, to parameterize 
the effect at the detector we need not only a magnitude of an elementary process giving a lepton flavor violation but 
also a parton distribution function which gives the dependence on the neutrino energy. We need a precise knowledge 
about nuclei so it is very difficult and we ignored these effects. 

Then we calculate how many no-oscillation neutrinos are necessary to observe the effect of new physics considering 
only statistical fluctuation. For a neutrino factory we translated this neutrino number into the parent muon number. 
If we know all the mixing parameters exactly, then it will be very easy to observe an effect of new physics, which we 
parameterized as in eqs. and (|r^. The optimum setup for the baseline length and the energy region in a 
neutrino factory is easily understood by the high energy expansion of the transition rate as given in eq.(p^ and 
eq.(|^), for example. However, in general, we cannot expect that we know the mixing parameters exactly. If the new 
physics effect gives a similar energy dependence for the transition rate with that of the weak interaction as eq.(p^, 
we cannot expect that this effect is observable eventually. On the contrary, if the energy dependence is different as 
eq.(|^), we can expect, even if we have the uncertainties of the mixing parameters, the new physics effect makes a 
significant modification to the transition rate. In this case the uncertainties of the mixing parameters do not make 
the observability of the new physics effect worse. To survey the presence of e^’™, the use of the appearance channel 
of Va gives the best sensitivity. The sensitivities to new physics effect depend not only on their magnitudes but 

also their complex phases. More precisely we can observe an effect of a certain combination of the CP phase S and 
these parameters as given in for example eq.(^). The sensitivities are proportional to |ep. Naively we can expect 
their significant contribution as long as these combinations are larger than 0(10“^). 

Finally, we gave a discussion about possible new physics and their consequences to the parameters e®’™. Many 
models predict large lepton flavor violations to explain neutrino Majorana masses. These violations also affect decays 
of particles. For example, in many models there will be a flavor changing fj, decay. Also these effects will appear in 
the matter effects and possibly at the detection processes. Their effects must be taken into account in the analysis of 
oscillation experiments if we take those models seriously. Among these effects fx ^ t changing effect is expected to 
be large to explain the large mixing in the atmospheric neutrino anomaly. To observe these effect the disappearance 
channel of > Vt channel are very effective. This changing effect may play an important role in 

appearance experiments like OPERA. Q 

Incidentally we briefly mention the result of the LSND experiment, In general, its result is interpreted by the 
neutrino oscillation. If its result is partially confirmed by MiniBooNE 1^] , namely, MiniBooNE see a flavor changing 
signal with a small rate which is at the lower end of the LSND result, then we may interpret the result by the flavor 
changing decay of tt and/or If we take this interpretation, then the expected signal in a future neutrino-oscillation 
experiment is quite different from that by a four-generation model. We need to keep this possibility in mind. 
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1. up to 0 {Sm 2 i, e) 


2 2 -2 '^+ 
-I'm — *23*2x13 ' 


4K. 


1 


+ :rC5S2x23S2xl2S2xl3 






A_ +"'3*13-13;^ ,, 2 


Srun 


'21 


. 9 JmS, c ?9 — A_ ^ 

sm"’-- L 

4E,. 


Sm^i 


y_ + *13*13-13 ^2 


Sm^i 


sin 


+ CioC 


(5m|i 


Sm^i 


13^^13-13 ^^2^^2^ _ y_ *13*13-13 _ ^^2^^2 ' 


2 A+ Jw 21*^12 

4E, 

2 A+ — A_ 


4E„ 


1 / ^JTT'21 

+ ^SsS2x23S2x12S2x13 ^Cf3Cj3_f3 ^2 _ ^ *13*13-13' 


Sm^i 


X I sin + sin -^- L - sin 

2 £'iy aEu ‘2,Ejj 


(Al) 


+ 2S23S2X13 (c5R.e[e®^] + sMm[e®^]) 


9 , 9 ^ 171.91 Cl 9 A— 9.9 A-U Jm .91 Cl 9 

- £LJA - L - C23 sin'^ —-^A_t£ 


X I C23 sm 


4E.S 


+ S 23 S 2 X 13 (* 5 F!.e[e®^] - cMni[e®^]) ( c^g sin 


Smlj cl2 - A 


2A,. 


4A,. 


-L + C 23 sin 


P + *23'-2x13 

A+ - (5wiicf2 

2E.S 


2 X + — A_ 

4E,x 


L 


L “t“ ^23 


A+-A_ 

2E„ 


L 


- S23S2X23S2X13 (c5Re[e^^] + sMm[e®^]) (^sin^ 

- ^32332X2332X13 (*5R'e[e®^] - cMm[e^.,.]) ^sin 

,^33 f C5Re[ae^] - sMm[ae^] , A+ - A_ ^ 
+ 2*23*2X13 sm^^L 


^2 X fa 


sm 


(A2) 
2 A+ — A 


(5m9iC?9 — A_ ^ . A+ — SrriysCsn , , A+ — A_ 

--L + sm—t-£l—t£;^_sin ^ 


2E,s 


4E,x 


L 


(A3) 


+ S2 X 23 C 23 52x13 


2 C 5 Re[oe;i] - S5lm[ae^] 2 C 5 Re[ae^] - S5lm[a e/x 2 Sml^cl2 - A 

+3 Jm^c^-A- 13 1 X.2 1*111 


Aj- — Sm^s c? 


- 21^12 


4E,. 


f^2 CjRejQep] - 5^Ini[ae^] ^2 QRe[ae^] - SjIm[oe^] ^ ^.^2 A+ - Sml^cl 


"13 (5moiC?9 —A_ 


13 A+ — Snioi c?9 


sm 


‘2l'-12 


4E, 


+ c 


2 CjRe[ae;^] - s^Im[ae^] 2 QRe[ae^] - s, 5 lm[ae^] A . 2 A+ - A 

13 i5mniC?9 — A- 13 ' ■ — ^ 


A+ — Smni c ?9 


4A,. 


-A 


+ ^32x2302332x13 ( Gs 


2 5jRe[ae;j] + c^Im[oe^] _ 2 s^Re[ae;x] + C 6 lm[ae^] 


Jmficf 2 - A_ 


13 


A+ - (5miiC?2 














































21 


. (5m^ic?2 - A- ^ , . A+ - dml^cl^ ^ . A+ - A_ 

X 1 sm-- L + sin-- L — sin —- L 


2E, 


2E,. 


2E,. 


, ^ 3 2 f C5Re[ae^] - S5lm[ae;.] \ 2 A+ - A_ ^ 

+ 452352X13^2X13 I 4E, ^ 


+ 52352x235 


52352x2352x13 


3 / C5Re[aeT-] - S5lm[aex] A . A+ - A_ 

'2X13 Lj sin^^L 


^ aRe[aeT] - 5jlm[aer] 2 QRe[aeT] - 5aIm[aeT-] \ . 2 ^”^ 21^12 ~ A- ^ 

— X- A+— i5m23cf2 / 'iEi, 


' 21^12 

2 CjRe[ae-r] - 5aIni[Oer] 2 CsRelaer] - 5jlm[ae-r] \ . 2 A+ - SmliC ^2 T 

5moiC?9 —A_ A+ —^mniC? ~ 


+ c: 


13 \ „2 I 

C5Re[aeT] — S5lin[aeT-] 2 C5Re[aeT-] — 55lm[aer]\ . 2 A+ — A 


(5miTC?9 —A_ 4" 5^3 A+— (577123 c^2 


(- 2 li '12 
2 




S5R,G[Q-e'7-] H“ C^Iin[Q-e7-] 2 5(5R'^[^eT] H” C^Ini^CleT-j 


252352X2352X13 J^ 2 ^^ 2 ^ _ ^_ ^13 X+ - Snij^cl^ 


X sin • 


21^12 

(5777|iC?2 - A 


2E,, 


L + sin 


A+ - Smlj^cl2 r A+ - A_ ^ 

— 1 j — sin- J-j 


2E„ 


2E,. 


2 ( C5Re[aeT] - 55lm[aer] \ . 2 A+ - A_ 

+ 2s2352x23S2xf3C2xl3 ( -- ) Sm L 


2 2 ( RG[ 0 /iT] j \ . 

+ 52352x2352xf3C2xl3 ( /^Ey 


A4--A 


2E,, 


-L 


— 2'®2x23C2x2352xf3 


Re[( 


l/iTj 


Re[(: 


'/iTj 


Sml^cl2 - A_ A+ - 5ml^c\2 

Re[a3,r] ^ Re[a3.r] \ ._^2 A+ - Sm^j^cf2 


n Ct 2 A — 

sin"'-- L 


Sml^cl^ - A_ A+ - Sml^cl^ 

'^T 


+ 


'' 21^12 

Re[a 3 JT-] Re [a,. 


sin 


Sm2iCi2 — X- A+ — Sm2iC 


12 


4.E, 

Ai - A. 


4.E, 

L 


sin 


2 -^3- 


4E,. 


-L 


1 


Im[ 


Im[ 


^52x2352x13 ^^^2^^2^ — A_ A+ — (57Tl2iCf_ 

SfTini Cl n A— , A-i_ S'm^'x Cl 9 

-^ + ^-" 2 E, 


O 2 2 / R' 6 [a^T-j 

^S23S2x2352>^f3C2xl3 1 ^ ^ 


'//rj \ . 2 '^+ r 

^ ‘ sm -- L 


4E. 


2 2 
- 5oQ . . ^'■o C 


23*2xl3‘'2xl3 

R,G[flee 


Re[aee] sin + 2sW 


+ 2 sooS?..,-„c 


2X13^2X13 I A+-A_ 


sm 


2E, 

2 A+ — A_ 


, A+ — A_ 




23*2X13 I 4^ 


4E„ 


AE, 


+ 5lc.S^,,,9,C 


23*2x13‘'2x13 I 4^ 


sin + 2si,s^ 


1 2 2 

“ 4'®2x23*2x13 


Re [a 


'mmJ 


2E, 

Re[niiii] 


'“[“>‘>‘1x5 sin" ^ 


23*2x13 4 ^ 

2 ^^2l'"12 ~ A— 


(5mfiC?2 - A- A+ - (577 i|iCi2 


sm 


2 „2 ; 4 ^^ 


L 


AE,. 


L 


+ 


(577i|iC?2 - A_ A+ - 6ml^cl2 

Re[ni(p] Re[(imi] 


. 2 A+ - 5ml^cl2 ^ 

sm -—- L 


(577i2iCi2 — A_ A+ — 5m\^c 


sm 


4.E, 

2 x+ — X- 


4E„ 


(A4) 


(A5) 


(A 6 ) 


(A7) 




































































22 


L 2 2 

~Sr)^^r>S 


Im[c 


'MMJ 


Im[Q 


2x23^2x1-3 ^^2^ C?,-A_ ^ X^- 6 mlcl 

L + sin 


X sin 


(5m^iC?2 - ^- T , ■ ^+~ '^™ 2 iC ?2 


- 2 s 


4 „2 


2E, 

Re[c 


'MMJ 


23'S2xf3'^2xf3 ^ 

1 2 2 _ /R£[^ 

^*2x23*2xf3^2xl3 1 


Sin 


2E, 

2 A+ - A_ ^ 
A+-A 


. A+ — A_ 


L sin 


^2 2 

^'S2x23'®2xf3 


Re[aT 


2£;. 

Re[ai- 


L -\- — s ^ 


2 ■’2x23-2x1-3 I 4 ^ 


Ini[o 


5m|iC?2 - •^- - <^"J 2 lC ?2 


Sin 


2 ^^ 21^12 

4E^ 


L 


Re[aT 


Re[aT 


(5m|ic?2 - A- A+ - i5m|ic?2 


+ 


'' 21'-12 

Re[aTT] 


sin 


A+ - ^ 


2 2 ^+ 


Re[Q 


i5mnnc?9 —A- A+ —i5m9iC?r 


sin 


4K 

2 A + - A- ^ 
4^;. 


^2 2 


Iin[Q 


Iin[Q 


2x23-2x1-3 U„2 2^ ^ X+ - Sml 


c?9 — A- 


X I sin + sin - sin 


2E,, 


-L 2 2 

~ 2'®2x23'®2x1-3''2x1-3 


Re[arr] 
A+ - A- 


Aj- — SrUni C?9 


2E, 

■ 2 ~ T 

sin —-L, 


2E,, 


4E^ 


L sin 


A+-A 


2 ^+ 


4£;„ 


where A±, 0i3 are defined as follows: 
1 


A± = 


5m\i + 5 m\is \2 + « i V {(<5^31 ~ <^W23sf2)c2xi3 — a} + {Srn^i — 6 m 2 iSi 2 )^s 


-2x13 


tan 2013 = ®2xi3(^m|i - (5miiS?2) 


C2xi3(^m§i - (5m^isf2) - a' 


See Ref.@ for details of the calculation. 


(AS) 


(A9) 


(AlO) 

(All) 


2. AP^s^^^{eij,T} up to 0{5m2i, e) 


AP,, 


= -2Re[e®^]s2x23 


C2x23Sf3 1 sin 


(5miicf2 - A 


+ C 2 X 23 sin 


4A, 

2 A+ - (5m|icf2 


-L — sin 


Aj- — 6m^^ c ?9 


4A„ 


iL + sin^ 


A+-A 


4P„ 


4P^ 


L - + C 2 X 23 ) sin' 


A+-A 


2 x\+ 


4P^ 


L 


- Im[e®]52x23 ( 5^3 sin 


Sm^Tciri — X-^ 2 . A+— 5m|icf2 


+ 2 


Re[c 


'flT\ 


4A„ 


L ) S 23 S 2 X 23 


'- 2 1'"12 

2E, 


2 2 

C 99 S 


-L — c^g sin ■ 


2E, 


L 


23 - 1-3 


(3 C2X13) sin 


(5miiC?2 - ^- 


'^ 23 ''! 3 (3 + ^2x1-3) 


2E, 

Aj — (5mni C ?9 


-L 


2E„ 


■-L 


A+-A- 


+ 2s2x23 


2'S23®1-3'^1-3‘^2 X 1-3 ■ 2E 


o Re[a;xr] 2 / 2 _ 2 „2 \ , r, Re[«/xr] 2 2 2 

(5m2 C ?9 - A_ "=2X23513(023 ^ 235 ^- 3 ) +2 ^ . 2 C2x23Cf3523Sf3 


-L 


(A12) 





















































23 


+ 


Re[i 


2 I , 2 '^^ 21^12 T 

Bm ——— L 


M-rJ g2 g2 


- I 2 


Re[ap 


-C2x23SfciSl3Cfci + 2 


Sml^cl^ - A. 


Re[a^ 


Z ( 1 Z Z \ 

^13^23^13 ^ _ gjj^2 ^2 '^ 2 x 23 C ; j ~ 3 ( C 23 — 8230^^) 


+ 


R/G^tt^T-j 2 2 2 1 • 2 SiTI‘2'^C 


A _ A_®23«2x 13C23 pm 


iE,, 


'■-L 


Re[Q 


Re[a^ 


j 2220 - —j AAA 

'5m2,c22-A_'''x23Sf3S23Cf3 " ^ _ ^^2^^2p2x23Ci3S23Sil 


I R-e[a/jT] 4 2 1 • 2 '^+ r 

+ A+ - A_ ^23^2X13^2X13 } Sin L 


(A13) 


[1] Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28 (1962) 870. 

[2] Super-Kamiokande Collab., Y. Fukuda et al, Phys. Rev. Lett. 81 (1998) 1562; Phys. Lett. B433 (1998) 9; Phys. Lett. 
B436 (1998) 33; Phys. Rev. Lett. 82 (1999) 2644. 

Kamiokande Collab., K. S. Hirata et al., Phys. Lett. B205 (1988) 416; Phys. Lett. B280 (1992) 146. 

IMB-3 Collab., D. Casper et al., Phys. Rev. Lett. 66 (1991) 2561; Phys. Rev D46 (1992) 3720; Phys. Rev. Lett. 69 (1992) 
1010 . 


SOUDAN-2 Collab., W. A. Mann et al., Phys. Lett. B391 (1997) 491; Phys. Lett. B449 (1999) 137; Nucl. Phys. Proc. 
Suppl. 91 (2000) 134. 

MACRO Collab., B. Barish et al., Phys. Lett. B478 12 0001 5: Nucl. Phv s. Proc. Suppl. 91 (2001) 141. 

N. Fornengo, M. Maltoni, R. Tomas, and J.W .F. Valle, tiep-ph/0108043 . 

[3] K2K Collab., J. E. Hill et al., ^iep-ex/0110034 , 

[4] MINOS Collab., P. Spentzouris et al , Nucl. Phys. Pr oc. Suppl. 100 (2001) 197. 

ICARUS Collab., F. Arncodo, et al., 5ep-ex/0103008. 

[5] OPERA Collab., A. Rubbia et al., Nucl. Phys. Proc. Suppl. 91 (2000) 223. 

[6] SNO Collab., Q. R. Ahmad et al., Phys. Rev. Lett. 87 (2001) 071301. 

Super-Kamiokande Collab. Y. Fukuda et al., Phys. Rev. Lett. 82 (1999) 1810; Phys. Rev. Lett. 82 (1999) 2430. 

GNO Collab., N. Ferrari et al., Nucl. Phys. Proc. Suppl. 100 (2001) 48. 

GALLEX Collab., W. Hampel et al., Phys. Lett. B447 (1999) 127. 

SAGE Collab., J. N. Abdurashitov et al., Phys. Rev. C60 (1999) 055801. Nucl. Phys. Proc. Suppl. (2001) 91 36. 
Kamiokande Collab., Y. Suzuki et al., Nucl. Phys. Proc. Suppl. B38 (1995) 54. 

Homestake Collab., B. T. Cleveland et al., Nucl. Phys. Pro c. Suppl. 77 (199 9) 13. 

J. N. Bahcall, M.C. Gonzalez-Garcia, an d C. Pena-Gara v. hep-ph/0111150; G. L. Fogli, E. Lisi, A. Palazzo, F. L. and 
Villante, Phys. Rev. D63 (2001) 1130 16 [tiep-ph/0102288|; G.L. Fogli , E. Lisi, D. Mon tanino, and A. Palazzo, Phys. Rev. 
D64 (2001) 093007 |hep-ph/0106247|; M. V. Garzelli, and C. Giunti, hep-ph/0111254. 

[7] L. Wolfenstein, Phys. Rev. D17 (1978) 2369; S. P. Mikheev and A. Yu. Smirnov, Sov. J. Nucl. Phys. 42 (1985) 913; H. 
A. Bethe, Phys. Rev. Lett. 56 (1986) 1305; V. Barger, K. Whisnant, S. Pakvasa, and R. J. N. Phillips, Phys. Rev. D22 
(1980) 2718. 

[8] KamLand Collab., A. Piepke et al., Nucl. Phys. Proc. Suppl. 91 (2001) 99. 

[9] Borexino Collab., E. Meroni et al., Nucl. Phys. Proc. Suppl. 100 (2001) 42. 

[10] LSND Collab., C. Athanassopoulos et al., Phys. Rev. Lett. 81 (1998) 1774; A. Aguilar, et al., hep-ex/0104049; K. Eitel, 
New J. Phys. 2 (2000) 1.1. 

[11] MiniBooNE Collab., A. Bazarko et al., Nucl. Phys. Proc. Suppl. 91 (2000) 210 hep-ex/ 000905f| . 

[12] Chooz Collab., L. Mikaelyan et al., Nucl. Phys. Proc. Suppl. 91 (2001) 120; Phys. Lett. B466 (1999) 415. 

[13] Palo Verde Collab., F. Boehm et al., Phys. Rev. D64 (2001) 1 12001. _ 

[14] H. Minakata, and H. Nunokawa, Phys. Lett. B504 (2001) 301 hep-ph/0010240| . 

[15] J. J. Go mez-Candenas, A . Blondel, J. Burguet-Castell, D. Casper, M. Donega, S. Gilardoni, P. Hernandez, and M. 


Mezzetto, tiep-ph/0105297. 

V. Barger, D. Marfatia, and K. Whisnant, hep-ph/011211£. V. Barger, ef al., hep-ph/0103052 


Raja, and K. Whisnant. Phv s. Rev. D63 (2001) 113011 |hep-ph/0012017 ]. 


V. Barger, S. Geer, R. 


B. Richter, hep-ph/0008222 


Phy s. Rev. D55 (1997) 322 lhep-ph/9605413ll. 


M. Tanimoto, Phys. Lett. B435 (19981 373 I hep-pli/98063TO 
H. Minakata, and H. N unokawa, hep-ph/ 0111130; hep-ph/0111131; Phys. Lett. B495 (2000) 369 |hep-ph/0604114 
Rev. D57 (1998) 4403 lhep-ph/97052 08l: Phvs. Lett. B413 (1997) 369 |hep-ph/970628l|. 

JHF-Kamioka project, Y. Itow et al., ' piep-ex/0106019| . 


Phys. 






























































24 


[16] 


N. Okamura, M. Aoki, K. Hagi wara. Y. Havato, T. Kobayashi, 
Study G roup. Y.F. Wang et. al, tiep-ph/0 11 13171 . 


T. Nakaya, and K. Nishikawa, [iep-ph/010422C; VLBL 


J. Sato, hep-ph/0008056 
119905 


J. Arafune, M. Koike, and J. Sato, Phys. Rev. D56 ('199 7') 3093: erraturn ibid. D60 (1999) 


hep-ph/97072TO 


J. Arafune, and J. Sato, Phys. Rev. D55 (1997) 1653; M. Koike, and J. Sato, , _ 

Neutrino Fa ctory and Muon Collider Collab., D. Ayres et al., physics/9911009; C. Albright et al., hep-ex/0008064; fiep- 
ph/0111030| . 


V. Barg er. S. Geer. R. R aia. and K. Whisnant, Phys. Re v. D63 (2001') 03 3002 hep-ex/000602f j; Phys. Rev. D62 (2000') 
013004 ||hep-ph/9911524|| ; Phys. Rev. D62 (2000) 073002 ||hep-ph/0003184|l; Phys. Lett. B 485T^00) 379 ||hep-ph/00042o4 
; V. B arger. S. Geer, a nd K. 'Whisnant, Phys. Rev. D61 (2000) 053004 |hep-ph/9906487 ; S. Geer, Phys. Rev. D57 (1998) 
6989 |hep-ph/971229C|, erratum ibid. D59 (1999) 039903. J. Burguet-Cast ell. M. B. Gavela, J. J. Gomez-Candenas, P. 
Hernandez, and O. Mena, Nucl. Phys. B608 (2001) 301 hep-ph/0103258| ; A. Cervera. A. Donini. M. B . Gavela, J. J. 
Gomez-Cadenas, P. Hernandez, O. Mena, and S. Rigolin, Nucl. Phys. B579 (2000) 17 |hep-ph/000210f j; erratum ibid . 
B593 (2001) 731; A. Donini, M. B. Gavela, P. Hernandez, and S. Rigolin, Nu cl. Phys. B^4 ( 20001 23 lhep-ph/9909254l: 
A. De Riijula, M. B. Gavela, and P. Hernandez, Nucl. Phys. B547 1199 9') 21 lhep-ph/981 139C]. 

M. Freund. P. H uber, and M. Lindner, Nucl. Phys. B615 (2001) 331 ||hep-ph/010507l|| ; Nucl. Phys. B585 (2000) 10 5 
[hep-ph/0004085|; M. Freund, M. Lindner, S. T. Petcov, and A. Romanino, Nucl. P hys. B578 120001 27 hep-ph/9912457|; 
K. Dick, M. Freund, M. Lindner, and A. Romanino, Nucl. Phys. B562 (1999) 29 |hep-ph/990330f j; A. Romanino, Nucl. 


Phv s. B574 12000') 6 75 |hep-ph/990942E|; E. K. Akhmedov, P. Huber, M. Lindner, T. Ohlsson, Nucl. Phys. B608 (2001) 
394 ||hep-ph/0105029| . 

M. Campa nelli. A. Bueno, a nd A. Rubbia, |hep-ph /9905240; A. Bueno, M. Campanelli. a nd A. Rubbia, Nucl. Phys. B589 


(2000) 577 |hep-ph/0005007|; Nucl. Phys. B573 (2000) 27; A. Rubbia, hep-ph/0106088 
P. Lipari, Phys. Rev. D64 (2001) 033002 fhep-ph/0102046i. 


M. Tanimoto, Phys. Lett. B462 (1999) 11^ hep-ph/9906516 
H. Yokomakura, K. Kimura, and A. Takamura, Phys. Lett. B496 (2000) 1’7 5 [|hep-ph/00091B 


J. Pinney, and O. Yasuda, Phys. Rev. D64 (2001) 093008_[^iem£li^[01050OT |; O. Yasuda, |hep-ph/0005134 

. “ ~ jP 


[17] 

[18] 


[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 
[23] 


[24] 

[25] 

[26] 

[27] 

[28] 

[29] 

[30] 

[31] 

[32] 

[33] 

[34] 

[35] 


[36] 

[37] 

[38] 

[39] 


M. Koike, and J.Sato, Phys. Rev. D62 (2000) 073006 [|hep-ph/9911258|; Phys. Rev. D61 (2000) 073012 ||h ep-ph/990946S|[ 
erratum ibid. D6 2 (2000) 079903; T. Ota, J. Sato, arid Y. Kuno, Phys Lett B520 289 |hep-ph/0107007 |; M. Koike, T. 
Ota, and J. Sato, hep-ph/0011387 to be published on Phys. Rev. D. 

M.C. Gonzalez-Garcia. Y. Gro ssman. A. Gusso. and Y. Nir, Phys. Rev. D64 (2001) 096006 Jiep-ph/0105159|; Y. Grossman, 
Phys. Lett. B359 (1995) 141 ||hep-ph/9507344i . 

A. M. Gago, M. M. Guzzo, H. Nunokawa, W. J. C. Teves. a nd R. Zukanovich Funchal Phys. Rev. D64 (2001) 073003 
[hep-ph/0105196[; J. W'. Valle, and P. Huber, hep-ph/010819S. J. W. Valle, Phys. Lett. B199 (1987) 432; E. Roulet, Phys. 
Rev. D44 (1991) 935; M. M. Guzzo, A. Masiero, and S. T. Petcov, Phys. Lett. B260 (1991) 154; L. Wolfenstein, Phys. 
Rev. D20 (1979) 2634. 

S. Bergmann, Y. Grossman, and E. Nardi, Phys. Rev. D60 (1999) 093008 [hep-ph/9903517[. 

M. Koike, T. Ota, and J. Sato, in Ref.. 

T. Ota, and J. Sato, Phys. Rev. D63 (2001) 093004 I hep-ph /OOl 1234 1. _ _ 

A. Bueno, M. Campanelli, S. Navas-Concha, and A. Rubbia, hep-ph/oil2297; A. Rubbia, ^iep-ph/0106088 . 

M. Gell-Mann, P. Ramond, and Slansky, in Supergravity, ed. by P. van Niewenhuizen, and D. Z. Freedman (North-Holland, 
Amsterdam, 1979). 

T. Yanagida, in Proceedings of the Workshop on Unified Theory and Baryon Number in the Universe, ed. by O. Sawada, 
and A. Sugamoto (KEK, Tsukuba, 1979). 

R. N. Mohapatra, and G. Senjanovic, Phys. Rev. Lett. 44 (1980) 912. 

A. Zee, Phys. Lett. B93 (1980) 389. 

F. Borzumati, and A. Masiero, Phys. Rev. Lett. 57 (1986) 961. 

A. Ghosal, Y. Koide, and H. Fusaoka, Phys. Rev. D64 (2001) 053012. 

Y. Smirnov, and M. Tanimoto, Phys. Rev. D55 (1997) 1665. 

G. MaLaughlin, and J. N. Ng, Phys. Lett. B455 (1999) 224: E . Mitsuda, and K. Sasaki, Phys. Lett. B516 (2001) 47. 
Particle Data Grou 
LSND experiment 
here. 


D. E. Groom et. al, http://pdg.lbl.gov/; Euro. Phys. J. C15 (2000) 1. 


also shows a positive signal for neutrino oscillation which will be tested soon. 0 However we ignore 


Investigation about supernova may already have given the information about its sign.[[14| 

If the contribution from To or Wo is not dominant in eq.(^, we have to sum up appropriately. 

Eventually e// makes a shift on an observed rie experimentally. 

U" is not necessarily unitary. 

Strictly speaking, this quantity is “power of test” to distinguish two theories, the theory including new physics and the 
standard model. Furthermore, note that since we are interested in “goodness of fit” for two theories, we should compare 
this quantity with distribution function with one degree of freedom. More detail on statistics is found in Ref. [^ . 

We understand that systematic errors should be also taken into account. However in this paper we refer to the errors 
originated from statistics and we treat the errors in the theoretical parameters. 

To clarify the statement here we adopt large values for eejT*- The value of £e/i in Fig.^ and ^ are too large to avoid the 
current experimental bound. As we said first, are bound so strictly that these terms can not be observed. 

Including systematic errors we will have much benefit. 

This type of an interaction is also induced by a tensor particle exchange. 
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[40] In the authors use the very severe constraint for the nniversality. 

[41] Some of the bound indicated in Ref, fl^] are originated from taking into account the interference of the processes that 

the final states are different from each other, but such processes can not be added up quantum mechanically. We select 
the bounds given by the processes including the same final state. To derive the bound for we use the constraint from 
T fj,ee in Ref. p£|] . 


